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Abstract 

Using the Wick-Cutkosky model and an extended version (massive 
exchange) of it, we have calculated the bound states in a quantum field 
theoretical approach. In the light-front formalism we have calculated 
the bound-state mass spectrum and wave functions. Using the Terent'ev 
transformation we can write down an approximation for the angular de- 
pendence of the wave function. After calculating the bound-state spec- 
tra we characterized all states found. Similarly, we have calculated the 
bound-state spectrum and wave functions in the instant-form formalism. 
We compare the spectra found in both forms of dynamics in the ladder ap- 
proximation and show that in both forms of dynamics the 0(4) symmetry 
is broken. 

1 Introduction 

Dirac's paper on forms of relativistic dynamics Q made it clear that the real 
difficulty of constructing a Hamiltonian theory of interacting particles that sat- 
isfies the requirements of special relativity, is finding the proper form of the 
interactions. This problem can be solved in a natural way by resorting to co- 
variant field theory. Then the construction of all generators of the Poincare 
group can be found in any textbook on quantum field theory. However, the 
results given have usually a formal meaning only and much work needs to be 
done to turn them into useful formulas. In particular, drastic approximations 
need to be taken, which may or may not be justified. An approximation that 
at first sight looks very attractive is to expand the states of the sytem under 
consideration into Fock states and truncate the expansion at a reasonable point. 
Doing so, one obtains a quantum-mechanical many-body problem and can use 
the powerful! machinery that has been developed for many-body systems. 



The success of Fock-space methods depends crucially on the properties of the 
vacuum. Only if one can build on a simple vacuum the Fock-space expansion 
is useful. The condition that the vacuum be simple limits the choice of forms 
of dynamics in the framework of quantum field theory essentially to the front 
form, also known as light-front dynamics (LFD). In this form three components 
of the four momentum, p^, p^, and — (jP + p^)/y/2, are independent of 
the interaction, while p^ — {jp —p^)/y/2 contains interaction and is for this 
reason said to be a dynamical operator. The variable x'^ = (x" +x'^)/ ^1 that is 
conjugate to p", is the evolution parameter of states and consequently denoted 
as light-front time. The dispersion relation of energy and momentum for a 
particle with mass m takes the formp" = + w?) / {2p'^) , p± = {p^,p^). One 
sees immediately that states of positive and negative energy can be separated 
kinematically, as positive energy and positive p^ are strictly disconnected. This 
property is called the spectrum property. In instant-form dynamics (IFD) where 

is the dynamical component of the four momentum, the energy may be 
either positive or negative, independent of the three momentum. An immediate 
consequence of the spectrum property is that massive particles cannot be created 
from the vacuum in LFD. This is strictly speaking not enough to reduce the 
true vacuum to the perturbative Fock vacuum in LFD and usually, in order to 
make progress, one makes the additional assumption that zero modes, states 
where all particles have p+ = 0, are decoupled. In this paper we shall also make 
this assumption. For an extensive review of many aspects of LFD we refer to 
Brodsky et al. (2). 

Adopting LFD we still may consider other forms of dynamics, in particular 
IFD, for comparison. Even though it is difficult to justify IFD in a genuine 
field theoretical setting, to be distinguished from a situation where binding 
energies are very much smaller than particle masses, it remains to be seen what 
the quantitative differences are between the results obtained in LFD and IFD 
respectively. On the level of calculating invariant amplitudes in a Hamiltonian 
approach the situation is clear. One can obtain the time-ordered amplitudes, i.e., 
those that one calculates in LFD or IFD, by integrating Feynman amplitudes 
over the relative energy variable, k~ in LFD and in IFD. Schoonderwoerd et 
al. have shown in Yukawa theory that the difference between the covariant 
box diagram and the ladder approximation to it in LFD is much smaller that in 
IFD. As this result is obtained for the box diagram with external particles on 
mass shell and initial and final state on the energy shell, one cannot immediately 
conclude that in a bound-state calculation LFD in the ladder approximation 
would be closer to a covariant calculation than IFD. In this paper we show our 
results in the simplest possible case, the Wick-Cutkosky model (WC) ||]. (This 
model has been reviewed extensively by Nakanishi Q.) 

The WC-model is concerned with two scalar fields, $ and </>, of masses m 
and ji respectively, with a coupling g<^*<b(j). The original WC-model has /i = 0. 
We study the bound-state spectrum and wave functions for a range of coupling- 
constant values g. Because in LFD, contrary to IFD, the components and 
Ly of the orbital angular momentum operator L contain interactions, manifest 
rotational invariance is broken in LFD. This means that if a truncation in Fock 
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space is made, the observables, in particular the masses of the bound states, 
will depend on the orientation of the light front. Below we explain how this 
violation of rotational invariance manifests itself in our calculations. 

It is well known that the WC-model exhibits 0(4) symmetry, like the non- 
relativistic hydrogen atom. This symmetry shows up as a degeneracy of the 
levels of different orbital angular momentum I. We study the breaking of this 
symmetry due to the ladder approximation, which amounts to the truncation to 
the two- and three-body sectors in Fock space. The violation of 0(4) symmetry 
cannot be expected to occur in LFD only and we show here, for the first time, 
that it also occurs in IFD. 

The WC-model is popular because it avoids the complications of spin and 
an exact solution of the Bethe-Salpeter equation is known in the case that the 
mass of the two-body bound state vanishes. Moreover, Cutkosky Q gives nu- 
merical results for the spectrum in a range of values of the coupling constant. 
Recently, Mangin-Brinet et al. |^ calculated the ground state of the WC-model 
in so-called covariant LFD, also making the ladder approximation. Ji and Furn- 
stahl estimated the masses of the IS and 2P states in the WC-model using 
a variational Ansatz for the two-body wave function in LFD, for small values 
of the coupling constant. Ding and Darewych jsj also used variational tech- 
niques to calculate the ground state of the WC-model. In the end, they solve an 
equation with a kernel that does not depend on the mass of state, which may 
explain the fact that there masses lie much below the masses found either in 
the Bethe-Salpeter formalism or in the other Hamiltonian calculations. 

Several authors have extended the WC-model by adopting a nonvanishing 
mass /i for the exchanged particle. We mention in particular Cooke et al. 
Nieuwenhuis and Tjon and again Refs. [|[ ^. In all these works the ratio 
/z/m = 0.15 close to the ratio of the pion to the proton mass was taken, making 
these investigations relevant for the deuteron. 

We consider our work as a necessary step to the treatment of more realistic 
models, e.g., the Yukawa model. In this domain Glazek et al. [0 have done 
important pioneering work in LFD. Fuda et al. fl^ constructed one-boson- 
exchange models for the nucleon-nucleon and the pion-nucleon interactions. In 
the latter case an LFD calculation was compared to an IFD one. Contrary to 
what we do, these authors fit the parameters in the two distinct forms of dynam- 
ics separately to the experimental data, while we keep the model parameters 
fixed and compare the spectra. 

This paper is organized as follows. In Sec.|| we give the details of the model 
we use and the expressions for the mass operators. Next we write down the 
effective two-body equations valid when Fock space is truncated to the two- 
and three-particle sectors. In addition we do not take self-energy terms into 
account, which would require a discussion of renormalization, which we want 
to avoid here. We do not give a detailed derivation of our equations, as e.g. 
Glazek et al. have done so in much detail for the Yukawa model and it is 
easy to adapt their methods to the case of scalar particles. In Sec.|^ we turn 
to the delicate question how to estimate the orbital angular momentum of a 
state in LFD, knowing that rotational invariance is broken. For the purpose 
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of characterizing the states we rely on Terent'ev's transformation which is 
known to be exact for states on the energy shell For bound states, which 
are by definition not on shell, this transformation can only give tentative results. 
The next section contains the details of our numerical methods and in Sec.|^ we 
give the masses and the wave functions we found. Finally we discuss our results 
and draw conclusions. 



Derivation of bound- state equations in the used 
model 



We use the Wick-Cutkosky model to describe the dynamics of two scalar par- 
ticles of opposite charge exchanging a neutral scalar particle. In the original 
model by Wick and Cutkosky [|j the exchanged particle is massless. It is pos- 
sible to extend this model to a version where particles with a certain mass are 
exchanged. In our model we have taken bosons of equal mass m and an ex- 
changed particle of mass /x, which may vanish. The Lagrangian for this system 
is given by: 



1^ 



— 02 - .g(/)$*$, 



(1) 



where g is the coupling constant, <f> the charged field of the bosons of mass m 
and (j) the field of the exchanged particle. 

In both forms of dynamics we have used a field theoretical method to de- 
rive the bound-state equation from this Lagrangian. In LFD there is also the 
possibility to make use of the so-called explicitly covariant light- front formalism 
(eg. Ref. ^). The bound-state equations in both methods are written in 
different coordinates, but it can be shown via a coordinate transformation that 
the equations are the same. 

Through the energy-momentum tensor it is possible to find an expression 
for the Hamiltonian. In LFD the Hamiltonian is given by: 



50$*$ 



(2) 



Here = {d^,d'^) and the integration element [d'^x] = dec d^a;^ 
In IFD the Hamiltonian is given by: 



H 



d^x 



a°$*a°$ -I- v$*v$ + m2$*$ 



(3) 



Here the integration element d^x is defined as d^x — dx^dx^dx^, as usual. 
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Following the standard procedures explained in Ref. |^ for LFD and Ref. |^ 
for IFD, we can expand the free fields (j> and <i> as follows 



/{d^/fc} (a(fc)e~*'=" + flt (fc)e''=") , (4) 
<S>{x) = /{d3fc}(6(fc)e-^'="+d^(fc)e*'="). (5) 



Here a annihilates the exchanged neutral scalar particle, b annihilates the charged 
boson and d annihilates the boson with opposite charge. 

The difference in the free field expansion in both forms of dynamics lies in 
the integration element and in the interpretation of the operators a, b and d. In 
LFD the integration element is 

In LFD we define the integration element in momentum space [d^fc] in a similar 
way as we did in coordinate space: [d^fc] = dfc+d^fc^. In IFD the integration 
element is given by 

{d3fc} = , (7) 

^ (27r)3/ V2-B(fc, m) ^ ^ 

Here d^k is the usual integration element in IFD and the energy is given by 



E{k, m) = yk^+rn? . 

Substituting the free field expansions into Eqs. or (^, we can find an 
expression for the Hamiltonian. Using normal ordering, we can write down the 
Hamiltonian, which can be split into two parts, a free part and a part containing 
the interaction. In Eqs. (^) and (||) the expressions are given in LFD: 



9 /"mS-iuS,/! 0{k+~k 



x(a(fc - k') [b^{k)b{k') + d^{k)d{k')] + a^{k - k') [b^{k')b{k) + d^{k')d{k)](p) 
The Hamiltonians in IFD are given in Eqs. 1^) and (0): 

Ho = [d^k[E{k,m){b''ik)b{k)+d\k)dik)) +E{k,fi)a\k)aik)] , (10) 



(27r)3/2 J y/8E{k, m)E{k', m)E{k - k' , ^i) 

X [a{k - k') {b\k)b{k') + d\k)d{k')) + a\k - k!) {b\k')b{k) + d\k')d{k))\l) 

In both forms of dynamics we can write down a Schrodinger like equation. In 
IFD we have the well-known equation H\'^) = £\'^) where £ is the total energy 
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of the system and \^) is the total wave function. In LFD the Schrodinger 
equation is given by 



2P+P--Pl)\^) = Af2|ijf). (12) 



Here P^ = Pf^^^ + P^^^i is the total Hamiltonian, M the total mass of the system 
and I^E') is the total wave function. 

Since it is not possible to solve a bound-state equation in full Fock space, 
we need to make a truncation. Here we follow the same procedure as Glazek et 
al. Taking only the two and three particle sectors into account, we write 
the total wave function as 

\^) = 1^2) + 1^3). (13) 
Using this truncation we can write the Schrodinger equation as 

H32 H33 J \ ip3 J \ ^3 J 

Here H22, H23, H32 and H33 connect a two particle state to a two particle state, 
a three particle state to a two particle state, etc. 

Using this Fock-space truncation and projecting out the states gives two 
equations, one for the two-particle sector and one for the three- particle sector 

{^2^22^2) + {V2\H23\^fz) = {ip2\E\ip2) , (15) 
{v^3\H32\(P2) + {f3\H33\^3) = {f3\E\f3)- (16) 

Eliminating either the two- or three-particle sector results in a bound-state equa- 
tion. We have chosen to eliminate the three-particle sector and write everything 
in terms of the two-particle wave function, which results in 

{E-H22)\^2) - H23 ^ ^„ H32W2). (17) 

The kets \(p2) and \(p3) in Eq. ( p^ can be written in terms of a wave func- 
tion and creation operators acting on the vacuum. We have chosen to use two 
different bosons in the two- and three-particle states. Then we find richer spec- 
tra, as symmetrization would remove anti-symmetric states. In the light-front 
formalism the wave functions are given by 



1^2) = j[d^p]^2{p,P-v)b\p)d^{P-pm. (18) 
|<^3) = [[d?p][d^p']^3{p,P-p',p' -p)b\p)d^{P-p')a\p' ~p)\0).{l9) 



Inserting these expressions into the bound-state equation and truncating after 
the second order in the coupling constant leads to an equation with four different 
contributions to the kernel. Two of these can be associated with the self energy 
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q .1-1 



-p^, 1-x -pj^ . 1-x' -p^. 1-x 

(b) 



Figure 1: First order time-ordered diagrams in the ladder approximation for 
one particle exchange. 



\j2 Pl+m' 



and are ignored in this paper. Only the terms which can be associated with 
the exchange of a particle between different constituents are taken into account. 
This means that we are working in the ladder approximation and the bound- 
state equation in this case becomes 

x{l-x)\ '''^^^-^'^^""(2^ jd-'^p'L^x'K{p±,x;p'i^,x')ip2{p'L,x') . (20) 

Here iir(pj_, a;') is the kernel of the equation, g the coupling constant 
and the relative variables x and p± are defined as: x — Pi /{Pi P± = 

(1 — x)pi^^_ — xp2^±. (Note that we used an imaginary coupling to obtain an 
attractive interaction.) 

The two terms present in the kernel of the bound-state equation can be 
represented graphically by two time-ordered diagrams (Fig. |l|). The two time- 
ordered diagrams arise due to the condition that the plus-momentum of the 
exchanged particle should be larger than zero (p'+ — p''" > or p+ — p'+ > 0). 
Each of these diagrams corresponds to an energy denominator which is present 
in the kernel of the bound-state equation. The full kernel is given by 

„ 1 f e{x'-x) e{x-x') \ ^, 

Here the theta-function gives the time-ordering and Da and Df, are the energy 
denominators corresponding to left hand side and right hand side graph in Fig. ^ 
The expressions for these energy denominators are: 

Da = pj+m^ p'^^n? {p'^-p\_f+p? ^22) 

X 1 — x' x' — X ' 

A = M.^P'.'^--' _pI^-^' _iP^~P'.?^^^\ (23) 
x' 1 — X X — x' 

In IFD the kets \ip2) and {ip^) in Eq. (|l3[) can be written as 



1^2) = Jd'p^2ip)bHp)dH~Pm, (24) 
1^3) = [d'pd'p'^3{p,p')bHp)dH-p'W{p' -pM- (25) 
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After doing calculations similar to those described above we get, when working 
in the centre of mass system, the bound-state equation in the IFD formalism 



9" 



(27r)3 



£ — ^ p '2 + m? — \/ p"^ + m? — {p' — p)^ + fj? 

^2{P') 



(26) 



4-\/p'^ + Ta^ \/ p"^ + \/ {p' — p)2 + /i^ 

The bound states in the ladder approximation are found by solving Eq. ( pO| 
and (EQ) respectively. 



3 Characterization of states 

After having found the bound states, we must identify the states and assign 
quantum numbers to them. This can be done by using the squared orbital 
angular momentum operator L ^ . In the IFD approach there is no problem to 
identify the states. The orbital angular momentum operator L is kinematical 
and both the orbital angular momentum quantum number / and the magnetic 
quantum number m are good quantum numbers. This in contrast to the LFD 
approach, where L is dynamical. This leads to the fact that only the helicity h 
is a good quantum number and that / is not. 

It is possible to approximate the L-operator in LFD by first transforming 
the light-cone variables (x,pj_) into the variables {pz,P±)- Note that this object 
is not a true three-dimensional vector, as its components do not transform 
properly under all 3D rotations. The transformation we use was first introduced 
by Terent'ev and is given by 



\/rri^ + W + Pz 



(27) 



Note that this transformation is exact for free particles. Its inverse is given by 



p1 



1/2 



1-x 



1/2 



.(28) 



After making this transformation, we can approximate the orbital angular mo- 
mentum operator L by using 



ip X Vp 



(29) 



To characterize a state we have found, we calculate the percentage of a specific 
angular momentum state (i.e. S, P, D, ...) present in this state. This is done 
by taking the inner product with a spherical harmonic Yjm, thus projecting 
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the wave function on a radial function. In this way we can calculate for every 
angular momentum quantum number I whether it is present in a calculated state 
and if so, how big its contribution to that state is. If the overlap is over 85% we 
characterize the calculated state with the quantum number We realize that 
this way of determining the angular momentum is only approximate. 



4 Method of solution 

We solved both bound-state equations (in LFD and IFD) numerically. This was 
done by integrating out the angular dependence and making an expansion into 
basis functions. 

The LFD wave function (p2{p±,x) in the bound-state equation, Eq. (po|), 
depends on both the momentum and the orientation of p± in space. It is possible 
to expand this wave function in cigcnfunctions of , 

(fi2{p±,x) = Y] nr-r-2 i^h{p±,x)xh{<P) ■ (30) 

V x(l~x) 



where h is the helicity and the factor 1/ y ~ -^'-^"^ introduced to sym- 

metrize the integral equation. The angular part of a wave function with fixed 
helicity h is 

Xhi(f>) = ^= exp(-i/i0) . (31) 

V ZTT 

We use it to integrate out the angular dependence in the integral equation, 
Eq. (|0|). 

As we can see from Eq. (|o|) the momentum dependent wave functions in 
LFD depend on both the perpendicular and plus component of the momenta. 
For the a;-dependcncc of this wave function we have used cubic spline functions 
Ref. |l^ and for the pj^-dependence we have used a basis which contains 
Jacobi polynomials, viz 

MP.) = ^MP^(^^T''pr''('4^)- (32) 

\P^ + KcJ \P +KcJ 

Here /i^c is a scaling parameter, I — 0,1,2,..., and Nki is the normalization 
constant, which is given by 



- 2V^,,fc!(fc + 2/ + l)! 
= T{k + l + l) ■ 

These functions are similar to functions which were first introduced by Olsson 
and Weniger (see below). We have chosen them because they have the right 
behaviour at the origin and for p —^ oo. 
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For technical reasons we symmetrized the integral equation using the trans- 
formation given in Eq. (^0[). The kernel of the symmetrized bound-state equa- 
tion for fixed helicity h is given by: 



2tt9{x' - x) 



3/2 /3/2 

Pi P± 



^jN'Px{l -x)- {p± + m?)^A'Px'{l - x') - + ^2) 

\h\ 



1 



+ 



2{x' - x)D'^ 



2tt0{x - x') 



1 



(x'-x)D'^ 



2{x - x')D'f^ 



( ■^p±p'^ V 

\(x-x')D'^ ) 



'^P±P'± 



)DL 



(34) 



where Z?^ and I?^ are given by 



D' = - 



p'l + rr? p'l + p\ 



1 - x' 



\-x 



(35) 
(36) 



Since we have integrated out the angular dependence the vector character of 
the momenta in these equations has disappeared. Note that in the kernel the 
absolute value of the helicity is present, which implies that states with opposite 
helicities are degenerate. 

The wave functions in IFD depend on the three-momentum (see e.g. Eq. (|26 
In this case we can expand the wave functions into spherical harmonics 



<^2(p) 



^Y^yim{p)(i>i{p) ■ 



(37) 



Note that we can integrate out both angles Q and (/>, whereas in LFD we could 
only integrate out one angle, 0. For the momentum dependent part of the wave 
function we can make an expansion into basis functions as well. Here we use a 
basis which was first introduced by Olsson and Weniger fia] 



Nkip' 



p^ + 



'+2 / 2 2 

(i+3/2,i + l/2) / P - ^J'sc 



pV 



p2 -j- 



(38) 



Here ^sc is a scaling parameter, Z = 0, 1, 2, and the normalization constant 
Nm is given by 



Nki 



2y/fisckl{k + 2l + 2) 
T{k + 1 + 3/2) 



(39) 
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We symmetrize the integral equation by using the transformation 



(40) 



After symmetrizing the equation and integrating out the angular dependence 
the kernel of the bound-state equation in IFD becomes 



Here I is the angular momentum quantum number and Vi{p^p') is obtained by 
integrating over the angle between p and p' and is given by 



where Pi{p ■ p') is a Legendre polynomial. 



5 Numerical results 

First we studied the accuracy of the matrix elements. Using enough integration 
points allows us to calculate the matrix elements with an accuracy of at least 
six decimal places. We studied the convergence of the spectra with respect to 
the number of basis functions, i.e., spline functions and Olsson Weniger like 
functions in LED and the functions introduced by Olsson and Weniger in the 
case of IFD. In LFD we found that 14 spline functions or more and eight Olsson 
Weniger like functions or more give an accuracy of at least three decimal places 
for the lowest state for any I. The states with principal quantum number n = 3 
have an estimated absolute error of 0.004. In IFD taking 20 basis functions or 
more gives an accuracy of four decimal places in all cases considered. 

Taking more basis functions into account gives some variation in the absolute 
values of the masses. The relative position of the bound states does not change, 
which means that the spectra are not changed qualitatively. 

We have also compared our calculations in LFD with those of Mangin-Brinet 
et al. 1^. They have calculated the ground-state masses only and we have 
compared our calculated ground-state masses with those given in Ref. We 
have to remark that in that paper a coupling constant a is used, where a = 
/ [IQttw?) . In Table [l] our calculated ground-state masses are compared with 
the values from Ref. |6| for different values of the coupling constant. For smaller 
values of g we find the same bound-state masses, while for larger values we find 





1 



(42) 



X 
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Table 1: Comparison of the ground-state masses given in Ref. (M.B.) and 
those obtained in the present calculation for different values of the coupling 
constant. 



a 


9 


mass M.B. 


mass this calc. 


6.0 


17.36 


0.820 


0.811 


5.0 


15.85 


1.160 


1.156 


4.0 


14.18 


1.415 


1.412 


2.0 


10.03 


1.790 


1.788 


1.0 


7.09 


1.923 


1.922 


0.5 


5.01 


1.973 


1.973 



somewhat smaller masses. As in Ref. [| the masses are given in three decimal 
places only, we conclude that our calculations are essentially in agreement with 
those of Mangin-Brinet et al. 

5.1 Spectra 

For different values of the coupling constant we have calculated the spectra in 
the case of massless exchange. In the LFD formalism we have calculated the 
bound states for different helicities and afterwards determined the corresponding 
orbital angular momentum quantum number. Here we have used eight basis 
functions for the p^-dependence and 14 spline functions for the plus-momentum 
dependence of the wave function. In the IFD formalism the bound states are 
calculated with / as a good quantum number and by using 25 Olsson Weniger 
functions. 

In Fig. 1^ the spectra are plotted for coupling constant g — 17.36, 14.18, and 
10.03, which correspond to a = 6, 4, and 2. We have only plotted the cases 
/ = 0, 1,2. We should remark that in Fig. || the spectrum for a = 6 is not 
completely shown. (We have left out the ground state to get a better view of 
the rest of the spectrum.) Due to the breaking of manifest rotational invariance 
in LFD the P-states and D-states are split. The center of gravity of these split 
levels is plotted in Fig. ^ as a dashed line and the values corresponding to this 
center of gravity are given in Table Looking at the S-states and the dashed 
lines in Fig. || we observe a more or less similar pattern as found in IFD (Fig. 
Looking at the spectra for different coupling constants we see that the binding 
becomes less when the coupling constant becomes smaller. This is the same in 
both IFD and LFD. Comparing between the different formalisms we see that 
the states in the LFD formalism have more binding. 

If instead of massless exchange 7^ is considered, one may expect the 
masses to increase. This is borne out by comparing fi = and /i = 0.15. The 
comparison is done for all three values of the coupling constant used above, i.e., 
a = 6, 4 and 2. In Table ^ the LFD calculated masses with orbital angular 
momentum quantum number / = 0, 1,2 (/i = and fi = 0.15) are compared 
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a = 6 
1 = 1 



a= 4 

1 = 1 



1 = 2 



a = 2 
1 = 1 



Figure 2: Spectra in LFD in the case = for different coupling constants. 



a= 6 

1=0 1=1 1=2 



a= 4 

1=0 1=1 1=2 



a = 2 

1=0 1=1 1=2 



Figure 3: Spectra in IFD in the case /U = for different couphng constants. 
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Table 2: Masses corresponding to the center of gravity for the P-states and 
D-states in LFD for fj, — and couphng constants g = 17.36, 14.18, and 10.03. 



g = 17.36 




g = 14.18 




g = 10.03 


I = 




I = 




/ = 


0.811 
1.744 
1.885 


IS 
2S 
3S 




1.412 
1.856 
1.936 


IS 
2S 
3S 




1.788 
1.945 
1.978 


IS 
2S 
3S 


I = 1 




I = 1 




/ = 1 


1.739 
1.883 


2P 
3P 




1.853 
1.934 


2P 
3P 




1.945 
1.977 


2P 
3P 


I = 2 




I ^ 2 




I = 2 


1.886 


3D 




1.935 


3D 




1.976 


3D 



to those calculated in IFD for a coupling constant g = 17.36 {a = 6). The 
calculated masses for a coupling constant g = 14.18 {a — 4) and g = 10.03 
(a — 2) are given in Table ^ and |^. In all these tables the quantum numbers of 
the states are given as well as the masses. 



Table 3: Comparison of the calculated mass spectra in LFD and IFD in the case 
of massless and massive exchange for a coupling constant g = 17.36 (a = 6). 













M = 


= 0.0 












LFD 


IFD 


LFD 


IFD 


LFD 


IFD 




I -- 


= 






/ = 


= 1 






/ = 


= 2 




0.811 


IS 


1.4931 


IS 


1.652 


2P 


1.8443 


2P 


1.858 


3D 


1.9259 


3D 


1.744 


2S 


1.8526 


2S 


1.782 


2P 






1.879 


3D 






1.885 


3S 


1.9304 


3S 


1.846 


3P 


1.9276 


3P 


1.907 


3D 














1.901 


3P 
























M = 


: 0.15 












LFD 


IFD 


LFD 


IFD 


LFD 


IFD 




I = 


= 






I 


= 1 








= 2 




0.923 


IS 


1.5384 


IS 


1.753 


2P 


1.9153 


2P 










1.846 


2S 


1.9276 


2S 


1.870 


2P 















As we already saw in Fig. the states have less binding for decreasing 
coupling constant and we find only a few bound states in the case oi g = 10.03. 
Tables |, |, and | show the same effect. From these results it also becomes 
clear that the states in the case of massless exchange have more binding than 
those for /i = 0.15. Especially in the case of Z = 1 or 2 the binding is sometimes 
so small that we did not succeed in calculating the bound-state mass with the 
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Table 4: Comparison of the calculated bound-state masses in LFD and IFD in 
the case of massless and massive exchange for a coupling constant g ~ 14.18 
{a = 4). 













^J■ -- 


= 0.0 












LFD 


IFD 


LFD 


IFD 


LFD 


IFD 




I = 


= 






I -- 


= 1 






I -- 


= 2 




1.412 


IS 


1.6666 


IS 


1.815 


2P 


1.9007 


2P 


1.923 


3D 


1.9534 


3D 


1.856 


2S 


1.9054 


2S 


1.872 


2P 






1.931 


3D 






1.944 


3S 


1.9558 


3S 


1.917 
1.943 


3P 
3P 


1.9543 


3P 


1.944 


3D 


















= 0.15 












LFD 


IFD 


LFD 


IFD 


LFD 


IFD 






= 






I 


= 1 






/ = 


= 2 




1.477 


IS 


1.7093 


IS 


1.898 


2P 


1.9628 


2P 










1.937 


2S 


1.9690 


2S 


1.945 


2P 















number of basis functions we have used. 
5.2 Wave functions 

Besides the spectra we calculated the wave functions corresponding to the bound 
states. In Fig. ^ the 2-D LFD wave functions of the IS, 2S and 2P states 
for g — 17.36 and ^ = are given. These figures show the variation of the 
wave function with x (or p+) and We see that the wave function goes 

asymptotically to zero when becomes large. 

In Fig. 1^ the wave functions corresponding to the IS, 2S, and 3S states 
are plotted for all three values of the coupling constant in the case of massless 
exchange. This figure shows that the wave function becomes flatter and broader 
when the coupling constant increases. All functions go asymptotically to zero 
and show the correct number of nodes; i.e. zero nodes for the IS state, one node 
for the 2S-state etc. 

In the case of the S-states, the helicity can only be equal to zero. This in 
contrast to the P- and D-states where we have helicity —1, and 1 (P-states) 
and —2, —1, 0, 1 and 2 (D-states). In our calculations we have worked with 
positive helicities only. We can do this because of the fact that only the absolute 
value of the helicity is present in the kernel, Eq. (|3^). In Fig. |6|the 2P, 3P, and 
3D states for all (positive) helicities are plotted in the case of /i = for coupling 
constant a = 6,4,2. As could be expected, the radial wave functions depend 
on the helicity, which is again due to the breaking of rotational invariance. The 
comparison between the radial wave functions calculated in LFD and IFD is 
made in Fig. |5.2|. In this figure the wave functions are plotted for a coupling 
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Table 5: Comparison of the calculated bound-state masses in LFD and IFD in 
the case of massless and massive exchange for a coupling constant g ~ 10.03 
{a = 2). 













^J■ -- 


= 0.0 












LFD 


IFD 


LFD 


IFD 


LFD 


IFD 




I = 


= 






I -- 


= 1 






I 


= 2 




1.788 


IS 


1.8474 


IS 


1.935 


2P 


1.9568 


2P 


1.973 


3D 


1.9801 


3D 


1.945 


2S 


1.9583 


2S 


1.949 


2P 






1.975 


3D 






1.978 


3S 


1.9808 


3S 


1.973 


3P 


1.9803 


3P 


1.978 


3D 














1.978 


3P 
























A* = 


= 0.15 












LFD 


IFD 


LFD 


IFD 


LFD 


IFD 






= 






I 


= 1 






/ = 


= 2 




1.833 


IS 


1.8849 


IS 


1.992 


2P 


1.9995 


2P 










1.995 


2S 


1.9982 


2S 


1.998 


2P 















constant g = 17.36 only. Comparing the plots in Fig. 5^ we see that the wave 
function calculated in the LFD formalism is flatter and broader than the wave 
function calculated in IFD. 



6 Discussion and conclusions 

We calculated the bound-state spectrum of the Wick-Cutkosky model in the 
ladder approximation both in light-front dynamics and in instant-form dynam- 
ics. We used as a variational method an expansion of the wave function in spline 
functions for the x dependence and an orthonormal set for the pj_ coordinate in 
LFD. Likewise, we used a similar set for the radial wave functions in IFD. We 
checked the absolute convergence of the spectra and also convinced ourselves 
that the relative positions of the mass eigenvalues will very probably not change 
at all when the number of basis functions would be increased above what we 
have used. 

As rotational invariance is broken in LFD by cutting off the Fock-space ex- 
pansion, we cannot expect the correct multiplet structure to appear in the ladder 
approximation. Moreover, because the angular-momentum operator is dynami- 
cal, we cannot immediately determine the angular momentum quantum number 
I of the states found in LFD. Using Terent'ev's transformation, which is approx- 
imate in the bound-state (off-energy-shell) case, we were able to characterize the 
bound states by calculating their overlaps with the spherical harmonics Y/m (p) . 
Following this procedure, we could estimate the degree of breaking of rotational 
invariance in LFD. Not surprisingly, rotational symmetry is less violated for 
small values of the coupling constant than for strong coupling. 
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Figure 4: LFD wave functions (IS, 2P, 2S) corresponding to the first three states 
in the case m = 1.0, ^ = 0.0 and g — 17.36 



The well-known 0(4) symmetry of the WC- model is violated in the ladder 
approximation in both LFD and IFD. If one computes the centers of gravity 
of the states with the same I and helicities h — one finds that the 

breaking of 0(4) symmetry in LFD and IFD are quite similar. 

In all cases we considered, the masses calculated in LFD are smaller than 
in IFD. One might try to explain this fact by pointing out that the Fock-space 
expansion converges more rapidly in LFD than in IFD, as shown in Ref. but 
one must keep in mind that in the latter paper the Yukawa model was investi- 
gated, so one cannot be absolutely certain of the validity of this argument for 
the WC-model. Calculations done by Cooke and Miller |^ and Schoonderwoerd 
et al. show that the breaking of rotational invariance in the Yukawa model 
is partially restored by including the box diagram. 

Although in the cases where the coupling constant is large, the bound sys- 
tems are certainly relativistic with binding energies of the same order as the 
masses of the constituents, we find even then that the binding energies follow 
the nonrelativistic rule En = Eijv? remarkably well. This feature is illustrated 
m Table |[ For very weak coupling, a = 0.5 and 1.0, we see a similar pattern 
for the Salpeter equation, see Table [7|, although in this case the binding is much 
larger than in the LFD case. 

The potential used is not bounded from below and shows a collapse of the 
wave function above a critical value of the coupling constant, see Ref. |19 . In the 
literature a critical value of a ~ 1.27 can be found Looking at Table we 
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Figure 5: LFD radial wave functions corresponding to the IS, 2S and 3S states 
for a = 6, 4 and 2. 

see that the ground-state mass for a = 2 has dropped enormously, illustrating 
the collapse. 

An analytical solution of the Wick equation in the weak-binding limit is 
given by Feldman et al. In Table || the bound-state masses are given for 

two values of the coupling constant, a = 0.5 and a = 1, calculated with the 
expression for the binding energy given in Ref . . For these two values of the 
coupling constant wc also calculated the bound-state masses in LFD and IFD 
(see Table ^) for Z = only. 

In his paper Cutkosky |^ gives a plot of his numerical results obtained in 
the WC-model. From his figure we read off masses for different values of the 
coupling constant (row ' WC in Table ^) . In a similar way we read off the masses 
from Fig. 2 in Ref. Q . Looking at this figure we see that only the results by Ji 
are not underestimating the WC-model masses. In Table ^ we give the masses 
which we read off from the curve labelled 'Ji'. For comparison we also give the 
masses of the IS state calculated in the LFD approach and those calculated 
with the analytical expression given by Ref. j2^ (row 'FFT' in Table H). Here 
we should remark that the analytical expression is only valid for small values of 
the coupling constant. 

From Table H we see that the masses calculated by Ji overestimate the bound- 
state masses for the WC-model. We also see that the analytical expression of |^ 
is only valid for small values of the coupling constant a. 

Besides the mass spectrum we have also calculated both the two-dimensional 
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Figure 6: LFD radial wave functions corresponding to the 2P, 3P and 3D states 
for a = 6, 4 and 2. 

and radial wavefunctions. Breaking of rotational invariance is clearly seen in 
the LFD case. 
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Table 7: Bound state-masses calculated with a relativistic Salpeter equation. 
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Figure 7: Comparison of the radial wave functions calculated in LFD and IFD 
at a coupling constant 5 = 17.36. 
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